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Wake potentials are derived for a closed (cylindrical) pillbox cavity as a sum over cavity modes. 
The resulting expression applies to on- and ofF-axis beams and test particles. The sum is evaluated 
numerically for a Gaussian drive bunch and compared to the wake potential derived from simulation. 



I. INTRODUCTION 

The wake potential describes the interaction between two charged particles mitigated by an external structure; 
the wake potential is especially useful for considering charged particles traveling at relativistic speeds on parallel paths 
through structures (such as radio-frequency accelerating cavities) in a particle accelerator. In this paper we derive 
the longitudinal wake potential for highly relativistic beams traveling parallel to the axis — but not necessarily on the 
axis — of a closed (cylindrical) pillbox cavity. 

The monopole wake potential for on-axis beams in a closed pillbox cavity was derived in Q (see also 0, 
Analytical wake potentials have also been found for other geometries, such as closed spherical cavities [1], conical 
cavities multi-cell pillbox cavities with elliptical cross section (for axial beams only) [6*1, and multi-cell pillbox 
cavities with beam tubes i8|] . We note that the wake potential for a pillbox cavity with infinitely long beam tubes 
differs significantly from the wake potential in a closed pillbox cavity: in a pillbox cavity with infinitely long beam 
tubes, the parts of the wake potential that vary azimuthally as cos{m9) vary radially as r™; in a closed pillbox cavity, 
the radial dependence is more complicated. 

II. WAKE POTENTIAL 

The wake potential describes the momentum change of a test particle caused by the fields excited by another charged 
particle. When a point charge qt {b for "beam") travels through a cavity, it creates an electric field E(r, 0, z,t) and 
also a magnetic field; the wake potential describes the momentum change of a test charge due to those fields. 

The wake potential is particularly useful when applied to highly relativistic particles, which maintain nearly constant 
speed while changing momentum. If a charge qt travels at light speed along the path (r, 9, z){t) — (r;,, 9b, ct), creating 
an electric field with z-component Ez{r,9, z,t), its longitudinal wake potential is: 

W,{s;rb,9b,rt,9t) = -— f dz E,{rt,9t,z,t=(z + s)/c). (1) 
qb J 

(We will consider a closed cavity of length £, so the integration runs from z = to z = ^.) If a test charge qt trails 
qii by a distance s, traveling along a parallel path {r,9,z){t) = {rt,9t,ct + s) from z = to z = its longitudinal 
momentum change due to the fields created by qb is 

Ap,{s; rt,9t) = -^W^is; r^, 9b, r*, 0*) (2) 
c 

(the test charge is assumed to be highly relativistic, so its speed remains constant even as its momentum changes). 
By convention, positive Wz corresponds to a loss in momentum. 

In cylindrically symmetric structures, it is helpful to decompose wakeflelds into azimuthal harmonics, i.e., fields 
with dependence cos{m9) and sin(m(?) for m = 0, 1, . . .: 

oo 

Wz{s;rb,rt,9t) = ^ Wz,.m{s;rb,rt) cos{m9t) (3) 

m=0 
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where we have chosen the beam to be at 9i, = 0; symmetry prohibits the excitation of any modes with sin(m6') 
dependence. In structures with cyhndrical symmetry and infinitely extended beam tubes, the azimuthal components 
of the wakefield have a particularly simple dependence on rt and (when they are smaller than the beam tube radius) 
i 

W,,^{s;n,rt)ocrTrr. (4) 

However, Ref. shows that this simple form depends on the fields at z = ±oo being the fields of a charge in an 
infinitely long beam tube; in a closed cavity, the fields at the ends of integration (at z = and z = £) upset the simple 
dependence on rj, and rt, though Wz^m still approximately assumes the above form for rj, and rt near the axis. 

III. PILLBOX WAKE POTENTIAL 

The monopole wake potential for a closed pillbox cavity has been derived for axial beam and test particles in [3]. 
In the same way, we derive an analytical expression for multipole wakes, for beam and test particles following paths 
parallel to the axis. We will write the wake potential as a sum over cavity modes; only TM modes (with magnetic 
field transverse to the z direction) will be excited, since the beam travels in the z direction and TE modes have 
Ez = 0. Choosing the beam to be at 6b = 0, only TM modes with cos(m6') dependence will be excited, and we need 
not consider modes with sin(m^?) dependence. 

We can write the field in the pillbox cavity of radius R and length £ as a. sum over modes; the (cosine) mode TMmnp 
(for integers m > 0, n > 1, p > 0) has fields (see [IJJ], Sec. 7.3.7): 

Ez,mnpir, 9, z, t) = COS (^) J„ cos(m^)e-*'^'"-* (5) 

Er,rnnpir, 9, z,t) = - (^) £- siu (^) (^^) cos(m0)e-^— * (6) 



Ee,mnp{r, 9, z, t) = sin [^^) J„ ) sin(m0)e-— ' (7) 

■2 

' sin(TO6l)e-''^""''* (8) 



TT ( a .\ ■ f^O^mnpR (p'KZ\m ( jm,nr\ 

Hr,mnp(r,9,Z,t) = I. COS I — 1 — J„ I I 

V Mo '^Jm,n \ I / T \ K / 

He,mnp{r,9,z,t) = ^ /lii^^ cos (^) j;„ cos(m0)e— (9) 

where Jm is the Bessel function of order m, J!^ is its derivative, and >n,n is the n zero of Jm- Mode TMmnp 
oscillates with frequency ojn 



mnp Jm^n ^ 



is 



The wake potential for a highly relativistic point charge traveling parallel to the z axis, at radius ri, (and 9b = 0), 

Wz{s-rb,ru9t) =2H{s)Y,K^np{rb,ru9t) COB {^^^^^) (11) 

vanp 

where H{s) is the Heaviside step function, and H{0) = 1/2, and kmnp{rb,rt,9t) is the loss factor for mode TMmnp- 

, , V*,„pirb,9b = 0)V,nnpirt,0t) 

kmnp{rb,rt,9t) ^ — (12) 

^^mnp 

where Kn„p(r, 0) is the (complex) voltage gain of a test particle crossing the cavity at transverse position {r,9), due 
to mode TMm„p when the cavity has stored energy Umnp in that mode: 

Vmnpir,9)= [ dzEir,9,z,t^ z/c) (13) 
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J„(:^ cos(™0)=^ (-l)''e- 



Jm.n 
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and 



1 + SmO ireoR^i ji 



mnp 



2-5po 4j 

where 6mn is the Kronecker delta. The loss factor is then: 



2 

7n,n 



2 mnp 



(14) 



2-5, 



'pO 



Jr. 



(^^^^ cos(m0t) • 2 [1 - (-l)Pcos(c^„„//c)] 



1 



^^O^jm.n Jm {jm,n) 



(15) 



The sum in Eq. pip unfortunately does not converge: for fixed m and n, the terms oscillate with constant amplitude 
as p increases. In case the beam and test particles travel alongthe same line, the sum can be analytically evaluated 
for s small enough that the walls at r = i? can have no effect [1, [![ — it sums to a sequence of delta functions, which 
offers some explanation for the sum's lack of convergence (considering the representation of a delta function as a 
non-convergent Fourier series). 

The sum's behavior can be improved by calculating the wake potential due to a charged bunch with total charge 
Qb and linear density profile A(s), rather than a point charge. The bunch wake potential is 



V4s)= / ds'Xis- s')W4s'). 



We will consider a Gaussian bunch 



A(.)^ 



27rcr 



exp 



s 

2^ 



(16) 



(17) 



and replace the cos{uJmnpS / c) term in Eq. (jlip . with the integral 



ds' 



■ exp 



1 



exp 



27rCT 



2a^ 



Re 



erfc 



cos(fcs') 

s + ia^k 
V2a 



(18) 



where k = comnp/c, and erfc is one minus the error function (of a complex argument). The bunch wake potential is: 



V.{s;n,rt,et) = 2iJ(s) ^ - 



2-d, 



pO 



m,n,p 



(^^) Jm (^^) cos{m0t) • 2 [1 - (-l)Pcos(w™„//c)] 



X - exp 
2 



mnp 



2c2 



Re 



V2ct 



(19) 



If we wish to know just the contribution from modes with cos(m0) dependence, we sum only over modes with that 



For a: ^ TO, Jm{x) « {x/2)"^ /ml |1^, so the multipole contributions to the wake are proportional to 



t 'b 



when 



r™ and r™ are small. 



IV. NUMERICAL PITFALLS 



The bunch potential in Eq. (|19p can be evaluated in a straightforward manner, except for the function 
g-y +2ja:iygj.fj,^_^ _ ^^^^ where we have used the abbreviations x = s/(\/2cr) and y = (jujranp/ {V2c). 

An easy and fast way to evaluate erfc uses the Faddeeva function w{z) = e^^ erfc(— iz) computed by the method 
of [H; i.e.. 



e 



-v^+^^^Verici-x - iy) = e""''u)(y - ix). (20) 
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However, w{y — ix) becomes enormous for ^ ?/^, while w{y — ix) remains tractable; similarly, for large y^, 
erfc(— X — iy) is enormous. In these cases we use the asymptotic expansion for erfc (see fl^l); because the expansion 
in for erfc(z) is valid for |argz| < 37r/4, we have to apply the identity erfc(— a; — iy) = 2~ erfc(a; + iy) before using 
the expansion when x > \y\. 

Specifically, to prevent overflow with double precision arithmetic, we perform a different evaluation for the following 
cases: if x^ ~ y^ > 500 and x > 0, we evaluate 



-v^+^i^Veiici-x - iy) w 26-^'+'^"^ 



(21) 



and if x^ — y^ > 500 and a; < 0, or if y^ — x^ > 500, then we evaluate 



-y-+2ixy 



erfc(— a; — iy) 



0F(x + iy) 



l-3...(2m-l) 

[2(a;2 - y2 + 2ixy)\' 



(22) 



These approximations yield nearly full accuracy in double precision when the specified conditions are satisfied. 



V. BEHAVIOR FOR LARGE n AND p 

In this section we show the behavior of terms in the sum in Eq. (I19|) for large n and p. For large p {ujmnp/c 3> 1/c) , 
we must consider two cases (the asymptotic limits are given in Sec. lIVP : the non-oscillatory contributions that depend 
on p behave as 



.p^/^erfc 



S + i(J^UJ,nnp/c 



\/2cr 

J 2e-'''('^'""''/'=)'/2cos(cj,„„ps/c) if s/cr > o-cj„„p/c 

When s ^ (7, convergence is very fast (with truncation error falling as the tail of a Gaussian), since the Gaussian 
charge distribution has had time to cancel out high-frequency contributions. When s is comparable to or smaller than 
a, the wake potential "feels" only part of the charge distribution, and higher frequencies matter more; consequently, 
the terms in the sum eventually fall off as l/a;^„p. 

For fixed m, but large n (n 3> mR/rf,, n 3> mR/rt), the envelope behavior (ignoring oscillatory contributions) of 
terms in the bunch potential series that depend on n is: 



3l,nJL{3m,n)^ (n-fm/2)%2^' 
decreasing slowly as ~ l/n?. 



^ (24) 



VI. WAKE POTENTIAL SIMULATION 



We compared the analytical wake potential against Wz,mis\ri,^rt) calculated via simulation, using the electromag- 
netic particle-in-cell (PIC) simulation capability of vorpal [14] in Cartesian coordinates. We excited a cavity using a 
current bunch at (and 0b = 0), with a Gaussian width a in the longitudinal direction z; the current bunch traveled 
at a highly relativistic speed and was (artificially) unaffected by the fields it generated. At each time step, we injected 
highly relativistic test particles (with charges too small to affect the cavity fields) at radius rt and regularly-spaced 
angles 0(, and recorded the momentum change of each test particle after crossing the cavity. We thus measured 
^Pz{s\rt^9t)^ which, after decomposition into azimuthal harmonics, yielded Wz^m(s; ''6, ?'*)■ 

With a radius R — 11.5 mm, the cavity's TMqio mode oscillated at 10 GHz; the length was chosen to be one-half 
wavelength at that frequency, £ — 15.0 mm. The Gaussian bunch [as in Eq. had a — 2£/25 (but was very thin 
in cross-section). Following the advice of Ref. [ll| (Sec. 3.2.3), we chose the cell length Az < y/a^Ji, resulting in 
Az = 0.333 mm and transverse cell sizes Aa:: = A?/ = 0.336 mm. 

The curved metal boundaries of the cavity were simulated using the Dey-Mittra algorithm [15l | , which requires (for 
improved accuracy) a reduction in time step from the standard Courant-Friedrichs-Lewy time step. 
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FIG. 1: Monopole wake potential (a) and dipole wake potential divided by the test-particle radius (b) in the closed pillbox 
cavity: the analytical (green line) and simulated (blue circles) wake potentials are nearly the same; the difference between 
them is shown by the magenta line; the red line at the bottom shows the drive current. The analytical sum was cutoff above 
Lojc — 300/(7 (/ ~ 1.2 X 10* GHz) for sja < 8.5 and lo/c — 10/a (400 GHz) for larger s. Simulated results are shown for cells 
of size 0.33 mm, and a time step half the Courant-Friedrichs-Lewy time step. 

Figure [1] shows the resulting m = and m — \ components of the wake (bunch) potential for identical beam and 
test-particle radii, — rt — 0.672mm. The difference between the analytical and simulated wakefield grows in time as 
the simulated modes slip in phase with respect to the analytical modes (due to error in the frequency of the simulated 
modes). 

Figure shows the wakefields compared to a higher resolution simulation, with w Ax = Ay = 0.17 mm. The 
error between simulation and theory is correspondingly reduced. 
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APPENDIX: FASTER CONVERGENCE WITH CROSS-SECTIONAL DISTRIBUTION 

The wakefield potential series can be made to converge faster (as n increases, and/or as m increases) by distributing 
the driving charge in the radial and azimuthal directions, as well as the longitudinal direction. 
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FIG. 2: Monopole wake potential (a) and dipole wake potential divided by the test-particle radius (b) in the closed pillbox 
cavity: the analytical (green line) and simulated (blue circles) wake potentials are nearly the same; the difference between 
them is shown by the magenta line. The analytical sum was cutoff above ujjc = 300/cr (/ 1.2 x lO" GHz) for s/cj < 8.5 and 
lij/c = 10/(7 (/ ~ 400 GHz) for larger s. Simulated results are shown for cells of size 0.17 mm, and a time step 40% of the 
Courant-Friedrichs-Lewy time step. 



To consider a distribution in 0, one has to add the beam angle d\, to the wake potential; this is done merely by 
replacing 6t ^ 9t ~ 0]y. One can then convolute the bunch potential with a normalized Gaussian distribution in 6^ to 
get an analytical result, as long as the distribution is narrow enough that the tails that extend beyond (0^ — tt, 0fc + tt) 
are negligible. 

To get an analytical result for a radial distribution is more difficult, but we have found a distribution that allows 
fairly convenient computation. We consider the function 



G(r;ro,(Tr) := 



exp 


ro(H- 


T-o/ffr) , 

r/ ar 

r 


exp 


-3^ 
r 


(r - ro)2' 




2rK 




+ ro/ar)ro/<7r^^ 


2rKo (2^(1 H 


h ro/ar)ro/(Tr j 



(A.l) 



{Km is a modified Bessel function of order m) with the following properties (for integrals, see Ref. 3.471.9 and 
6.635.3): 



Giro) = 



G'(ro) 
G"(ro) 



exp [-1 - 2ro /cr^)] 



2roKo (2^(1 + rn/crr)r()/cr, 



= 



1 + 2ro/cr, 



G{ro) 



(A.2) 

(A.3) 
(A.4) 



drG{r) = 1 



(A.5) 
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IG{k) 



drG{r)Jmikr) 



(A.6) 



IG{k) 
IG{k) 



^™ \te (1 + ^) 



J, 



K„ 



210. 



0-2 /!c2 



1 + 1 



a (^2^(1 + ro/(7r.)ro/o-r.^ 



(fej^o) when fc ^ l/ur and rp ^ ct^ 
. , s i^m ( v/2fc/cr,.ro ) 



■ exp 



2/<Tr — \/2kJor^ Tq 



when fc ^ 1 /(Tr and ro S> cr. 



(A.7) 



(A., 



That is, G{r; rg, tXr) peaks at r = tq, has unity integral, approaches 6(r — ro) for cr,. — > 0, and (most important) yields 
a not-too-complicated analytical result when integrated with Jm{kr). 

To find the bunch potential for a radial distribution given by G(r; tq, (7^) one merely needs to replace Jmijm.nfb/ R) 
in Eq. with Eq. (jA.6[) . where r — and fc = jm,nlR- 

There are a couple of numerical difficulties to consider. First, one should use the identity 



7 2 fj^ 
Jm.n r 



1-1 



■2 2 / d2 



'2 2 

J rn, 71^ r 



(A.9) 



i?2 



1 + 1 



in case jm,nO'r/R ^ 1- Second, Km(x) and Ko{y) become small enough to underflow numerical arithmetic as x and 
y become large, although their quotient may be of tractable magnitude. In this case one had better evaluate: 



Km{x) 



e^'Kmix) 



Ko{y) evKoiy) 
where the asymptotic expansion (see [T^]) is used to evaluate 



e^KAx) = ,l- 



1 + E 



(4t/2 - l)(4;y2 - 32) . . . (4^2 „ (2j ^ 1)2) 



(A.IO) 



(A.ll) 



for 1/ = 0, 1 (summing up to j = 5 yields nearly full accuracy for double precision when x > 400), and then recursion 
is used for higher ly; 



e^K,+x{x) = e"if,_i(x) - —e^K,{x) 

X 



(A.12) 
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